function [y, Position] = accuracyIndex(model, naive1, naive2, truevalue, R, C)
y = zeros(1, 3);
Position = zeros(R*C,2);
for i = 1:R
    for j = 1:C
        
        v(1) =  abs(model(i,j) - truevalue(i,j));
        v(2) =  abs(naive1(i,j) - truevalue(i,j));
        v(3) =  abs(naive2(i,j) - truevalue(i,j));
        
        [M,I] = min(v);
        switch I
            case 1 
                y(1) = y(1) +1;
                Position(i+(j-1)*7,1) = i;
                Position(i+(j-1)*7,2) = j;
            case 2
                y(2) = y(2) +1;
            case 3
                y(3) = y(3) +1;
        end
        
   end
end


end
